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On the morning of January 15, 2006, the Stardust capsule successfully landed at the Utah 
Test and Training range in northwest Utah returning cometary samples from the comet 
Wild-2. An overview of the entry, descent, and landing (EDL) trajectory analysis that was 
performed for targeting during the Stardust Mission Navigation Operations Phase upon fi- 
nal approach to Earth is described. In addition, how the predicted landing location and the 
resulting overall 99 percentile landing footprint ellipse obtained from a Monte Carlo analy- 
sis changed over the final days and hours prior to entry is also presented. The navigation 
and EDL operations effort accurately delivered the entry capsule to the desired landing site. 
The final landing location was 8.1 km from the target, which was well within the allowable 
landing area. 


I. Introduction 

S tardust, the forth of NASA’s Discovery class missions, was launched on February 7, 1999. The spacecraft per- 
formed a close flyby of the comet Wild-2 coming within 149 km of the comet nucleus. The cometary samples 
were collected by extending a collection tray on a boom into the gas/dust freestream emanating from the comet 
(Fig. 1), where the particles were trapped in a material 
called aerogel. Once the collection process was com- 
pleted, the collection tray was retracted back into the 
capsule. In addition to collecting cometary particles, 

Stardust also collected interstellar dust particles during 
its 7 year journey. Stardust is the first mission to return 
samples from a comet. Reference 1 gives an overview of 
the Stardust mission. 

Upon Earth return on the morning of January 15, 

2006, the entry capsule containing the cometary samples 
was released from the main spacecraft, and descended 
through the Earth’s atmosphere decelerating with the aid 
of a parachute for a successful landing at the U.S. Air 
Force’s Utah Test and Training Range (UTTR) in 
Northwest Utah. Reference 2 describes the Stardust 
Earth return trajectory strategy. Figure 1. Stardust spacecraft in sample configuration. 
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This paper provides an overview of the entry, descent, and landing (EDL) trajectory analysis that was performed 
for targeting the capsule to UTTR during the Stardust Mission Navigation Operations Phase upon final approach to 
Earth. In addition, how the predicted landing location and the resulting overall 99 percentile footprint ellipse (ob- 
tained from a Monte Carlo analysis) changed over the final days and hours prior to entry is also presented. This 
analysis was required in order to substantiate the robustness of the capsule descent to assure that all entry mission 
and public safety requirements were satisfied prior to gaining authorization for capsule separation from the main 
spacecraft for Earth entry. 


II. Capsule Overview 

The Stardust capsule (Fig. 2) is approximately 0.8 m in diameter. Its forebody is a blunted 60 deg half-angle 
sphere-cone. The afterbody is a 30 deg truncated cone. The entry velocity for the Stardust capsule was the highest 
(inertial velocity of 12.9 km/s) of any Earth returning mission to date. For comparison, the Apollo lunar missions 
had entry velocities of 1 1.0 km/s. This high entry velocity resulted in the highest heating rates for any Earth return- 
ing vehicle. Traditional carbon-phenolic based thermal protection systems (TPS) are very effective at such intense 
heating levels; however, they are quite heavy. To remain within project mass limits, a new lightweight heatshield 
material Phenolic Impregnated Carbon Ablator (PICA) was utilized. 2 3 Reference 4 provides an overview of the Star- 
dust capsule and its design. 



Figure 2. Stardust capsule configuration. 


III. EDL Overview 

Figure 3 shows the nominal Stardust entry sequence, with the terminal descent phases highlighted. Four hours 
prior to entry, the 45.8 kg capsule was spun-up to 13.5 rpm and separated from the main spacecraft. The capsule has 
no active guidance or control systems, so the spin-up is required to maintain its entry attitude (0 deg angle-of-attack) 
during coast. Throughout the atmospheric entry, the passive capsule relies solely on aerodynamic stability for per- 
forming a controlled descent through all aerodynamic flight regimes: free molecular, hypersonic-transitional, hyper- 
sonic-continuum, supersonic, transonic, and subsonic. The capsule must possess sufficient aerodynamic stability to 
overcome the gyroscopic (spin) stiffness in order to minimize any angle-of-attack excursions during the severe heat- 
ing environment. Additionally, this stability must persist through the transonic and subsonic regimes to maintain a 
controlled attitude at drogue and main parachute deployment. The inertial entry velocity and flight-path angle for 
Stardust were 12.9 km/s and -8.2 deg, respectively. 
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Spin rate = 13.5 rpm 



Figure 3. Nominal Stardust capsule entry sequence. 

Reference 5 provides an in-depth description on the development of the entry sequence; specifically, the use of 
the high spin rate and a supersonic drogue parachute. During descent, the capsule utilizes a g-switch (i.e., gravity- 
switch) and two timers for deployment of the drogue and main parachutes. The g-switch is triggered after sensing 
3 g’s (on the decelerating side). At that point, the drogue timer is initiated. After 15.04 seconds, the drogue para- 
chute is deployed (approximately Mach 1.37), and the main timer is initiated. After 350.6 seconds (approximately at 
an altitude of 3.1 km mean sea level [1.8 km above ground]), the main parachute is deployed to slow the capsule for 
landing at UTTR. This nominal entry sequence is sufficiently robust to accommodate off-nominal conditions during 
the descent as shown by the Monte Carlo dispersion analyses in Ref. 5. 

IV. Earth Return Strategy 

The Stardust event timeline for final Earth approach is shown in Figure 4, which highlights the Trajectory Cor- 
rection Maneuvers (TCM) that were baselined for attaining the proper entry conditions. Reference 2 provides an 
overview of the Earth return strategy showing all of the required TCMs. Prior to TCM- 18, which occurred at entry 
(E) minus 10 days (d), the Stardust return trajectory was on a path that missed the Earth. Only after TCM- 18 was 
successfully executed did the trajectory of the spacecraft become targeted to the Earth (within the atmosphere). 

Final targeting was accomplished with TCM- 19 at E-36 hours (hr), which placed the nominal landing location in 
the eastern portion of UTTR. If TCM- 19 had not executed or only partially executed, a contingency maneuver 
TCM- 19a or TCM- 19b would have been implemented at E-12 hours to achieve the desired target landing location. 
At E-4 hours, the capsule was separated from the main spacecraft, thus starting the EDL sequence illustrated in 
Fig. 3 The separation maneuver imparted to the capsule the remaining delta-v required to target the desired nominal 
landing location at the center of UTTR. At E-3.7 hours, a TCM was performed to divert the main spacecraft into an 
orbit ahead of the Earth. If TCM- 18, 19, or 19a or 19b had all been unsuccessful, the capsule/main spacecraft would 
have flown by the Earth as depicted in Fig. 4 by the solid line. During Mission Operations, both TCM-18 and TCM- 
19 executed very successfully, as did the separation and divert maneuvers. As a result, the desired entry conditions 
were achieved with amazing accuracy. 2,6,7 
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Figure 4. Stardust final Earth approach event timeline. 


V. Trajectory Simulation 
A. Entry Trajectory Requirements and Constraints 

The Stardust atmospheric entry trajectory was designed to fit within an envelope of derived requirements and 
physical constraints based upon the capsule hardware design. As such, for a successful landing, all entry requirements 
must be satisfied. Table 1 lists all the EDL requirements and their specific bounds. Monte Carlo dispersion analyses, 
described in subsequent sections, were performed during the Mission Operations Phase to assess the satisfaction of 
these requirements. 


Table 1: EDL Requirements and Constraints 


Requirement 

Limit 

Entry Flight-Path Angle Error, deg 

<±0.08 

Entry Attitude, deg 

< 10 

Max Heat Rate, W/cm 2 

< 1200 

Attitude at Max Heat Rate, deg 

< 10 

Max Heat Load, KJ/cm 2 

<32.0 

Max Deceleration, Earth g 

<40 

Drogue Chute Deployment Attitude, deg 

<30 

Drogue Chute Deployment Mach Number 

> 1.2 &< 1.6 

Landed Footprint, km 

<84 
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B. Monte Carlo Uncertainty Sources 

During the entry, off-nominal conditions may arise that affect the descent profile. These off-nominal conditions 
can originate from numerous sources: state vector uncertainties, capsule mass property measurement uncertainties; 
separation attitude and attitude rate uncertainties; limited knowledge of the entry day atmospheric properties (den- 
sity and winds); uncertainty with the aerodynamics; and uncertainties with parachute deployment. In the analysis, an 
attempt was made to conservatively quantify and model the degree of uncertainty in each mission parameter. For 
this entry, 41 potential uncertainties were identified. 5 Table 2 captures these uncertainty sources, along with their 
corresponding 3 -a variances. The subsequent subsections describe in greater detail a few of the uncertainty sources. 


Table 2: Monte Carlo Analysis Variables 


Variable 

3-a Variation 

Distribution 

Entry states 

Based on covariance (See Ref. 6) 

— 

Radial center-of-mass offset, mm 

±0.254 

Gaussian 

Axial center-of-mass, mm 

±0.254 

Gaussian 

Moments of Inertia (Ixx, Iyy, Izz) 

±2%, ±5%, ±5% 

Gaussian 

Cross products (Ixy, Ixz, Iyz), kg-m 2 

±0.003, ±0.003, ±0.003 

Gaussian 

Separation pitch and yaw attitude, deg 

±0.56, ±0.56 

Gaussian 

Separation pitch and yaw rates, deg/s 

±2.83, ±2.83 

Uniform 

Separation roll rate, deg/s 

±12.0 

Gaussian 

Aerodynamic coefficients 

See Ref. 5 

— 

Ablation mass loss 

±10% 

Gaussian 

Drogue Parachute Drag Coefficient 

±10% 

Uniform 

Main Parachute Drag Coefficient 

±10% 

Uniform 

G-switch acceleration trigger value 

±10% 

Uniform 

Drogue parachute timer, s 

±0.05 

Uniform 

Main parachute timer, s 

±0.05 

Uniform 

Atmosphere 

GRAM-95 model (See Ref. 8) 

— 


7. Entry Covariance 

The Stardust strategy for Earth approach was designed to maximize public safety in light of possible anomalies 
and contingencies, while still preserving the capability to meet the entry requirements. As a result, a series of ma- 
neuvers were performed to set up the approach and entry (see Ref. 2). Initial conditions at entry were obtained from 
orbit determination solutions performed by the Stardust Navigation Team. References 2, 6, and 7 provide a descrip- 
tion of the navigation process during the return phase, and the determination of the final arrival conditions prior to 
entry. The navigation accuracy obtained for Stardust yielded extremely small state errors upon Earth arrival. The 
final orbit determination solution produced an inertial entry flight-path angle of -8.21 deg with a 3 -a error of 
±0.0017 deg, which was well within the ±0.08 deg requirement. 

2. Capsule/Cruise-Stage Separation 

Based on the final main spacecraft and capsule mass properties, a statistical separation analysis was performed to 
predict separation attitude and attitude rate errors. The attitude errors predicted in pitch and yaw were ±0.56 deg in 
each axes. The attitude pitch and yaw rate errors were ±2.83 deg/s in each axes, and the roll rate error was ±12.0 
deg/s. These variations were used as inputs in the Monte Carlo analysis. 

3. Atmosphere Model 

The Earth atmosphere model utilized by Stardust for the entry trajectory design and analysis was the Global Ref- 
erence Atmospheric Model- 1995 Version (GRAM-95). 8 This model is an amalgam of three empirically based global 
data sets of the Earth that can produce an atmosphere profile as a function of altitude for a given date, time, and 
positional location about the Earth. GRAM-95 produces a representative atmosphere taking into account variations 
in diurnal, seasonal, and positional information for a given trajectory to produce nominal density, temperature, and 
pressure profiles and their statistical perturbations along the trajectory flight track. GRAM-95 is not a predictive 
model. A profile is generated based on historical data for a given time, season, and location. 
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Figure 5 shows samples of few randomly perturbed 
density profiles as a percentage of the nominal profile 
for the Stardust entry date of January 15, 2006 produced 
by the GRAM-95 model. Also, depicted are the upper 
and lower (±3 -a) boundaries of the possible density 
variation. As seen, density variations between ±5% to 
±25% are possible. In addition, GRAM-95 can also pro- 
duce nominal wind profiles and their statistical perturba- 
tions for the Northward, Eastward, and vertical wind 
components. Figures 6 and 7 show the nominal as well 
as a few randomly sampled wind profiles (for the Star- 
dust entry date of January 15, 2006) for the Northward 
and Eastward wind components, respectively, along with 
their upper and lower (±3 -a) boundaries. Also depicted 
is the actual wind profile obtained from balloon meas- 
urements 2 hours prior to entry (E-2 hr), which will be 
discussed in a subsequent section. In the Monte Carlo 
analysis, an atmosphere profile (density and wind com- 
ponents) was randomly generated for each case having 
the characteristics shown in Figs. 5-7. 

C. Trajectory Analysis 

Two trajectory propagation codes were utilized for 
the Stardust landing dispersion analyses: the Program to 
Optimize Simulated Trajectories (POST) program 9 and 
the Atmospheric-Entry Powered Landing (AEPL) pro- 
gram. 10 Two codes were employed to obtain independent 
verification that the predicted nominal landing location 
and the overall size of the dispersed landing footprint 
ellipse were within the UTTR boundaries to ensure pub- 
lic safety. 

The POST trajectory analysis was performed model- 
ing six-degree-of-freedom (6DOF) dynamics, which 
included all forces and torques on the spacecraft, from 
atmospheric interface to drogue parachute deployment. 
During this portion of the entry, the full set of capsule 
aerodynamics and mass properties were incorporated 
into the simulation to accurately model the hypersonic 
descent. 5 From drogue parachute deployment to landing, 
three-degree-of-freedom (3 DOF) dynamics were used, in 
which only the drag force was modeled and was as- 
sumed to act opposite the wind-relative velocity vector. 
The POST trajectory simulation seamlessly transitions 
from 6DOF to 3DOF dynamics within a single continu- 
ous simulation. 

The version of the AEPL program used for Stardust 
employed 3DOF analyses throughout. Since the Stardust 
entry was unguided and uncontrolled, the 3 DOF results 
from AEPL agreed well with the POST 6DOF/3DOF 
simulation. The POST results were baselined as prime 
for the mission. In general, there was very good agree- 
ment between the two simulations. 



U 1 1 — : 1 1 

0.7 0.8 0.9 1.0 1.1 1.2 1.3 

Density/Nominal Density 

Figure 5. Density variation from GRAM-95 model. 



Northward Wind Component, m/s 

Figure 6. Northward wind component variation from 
GRAM-95 model. 



Figure 7. Eastward wind component variation from 
GRAM-95 model. 
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D. Monte Carlo Dispersion Analysis 

A Monte Carlo dispersion analysis was utilized to statistically assess the robustness of the Stardust entry to off- 
nominal conditions to assure that all EDL requirements and constraints were satisfied (see Table 1). All the input 
variables listed in Table 2 were randomly varied in the Monte Carlo dispersion analysis within their respective vari- 
ances and distribution types. The analysis included uncertainties in the initial state vector, capsule mass properties 
(mass, center-of-gravity, inertia), separation attitude and attitude rates, aerodynamic coefficients, ablation mass loss, 
atmospheric density and winds, parachute drag, g-switch trigger value, and parachute deployment timers. 

For the dispersion analysis, 3000 random off-nominal cases were run for all the navigation orbit determination 
(OD) solutions that were computed 2,6 at the various event times during the Mission Operations Phase. This analysis 
was performed to determine the appropriate magnitude and direction for the TCM-18 and TCM-19 maneuvers for 
proper targeting to UTTR. In addition, this analysis was used to assess the OD solution stability and to understand 
the movement of the nominal landing location and the variation in the 99 percentile footprint ellipse size within 
UTTR. This understanding was crucial in order to gain authorization for capsule separation and the subsequent Earth 
entry. The size of the 99 percentile footprint ellipse obtained from the Monte Carlo dispersion analysis was used in a 
public safety probabilistic analysis to certify that the risks of the Stardust capsule entry were acceptable. Reference 
1 1 describes the hazard analysis that was performed for the Stardust entry capsule using the 99 percentile footprint 
ellipses generated by this Monte Carlo analysis. This hazard analysis was preformed for the nominal scenario of an 
intact capsule, as well as for a burn-up and breakup scenario in case of capsule failure during the entry. 

The Monte Carlo dispersion analysis was performed on all the post-TCM18 OD solutions (OD s06008a through 
OD s06014a). Figure 8 shows the corresponding results at landing assuming that a perfect TCM-19 is executed. For 
clarity, only the results for three OD solutions are shown (OD s06011d, OD s06012b2, and OD s06014a), where 
their nominal landing locations (center points) and the 99 percentile footprint ellipses at UTTR are depicted. The 
target location selected for Stardust is near the center of UTTR having the coordinates 246.55 deg East Fongitude 
and 40.3 167 deg North Fatitude. 



East Longitude, deg 

Figure 8. Landing locations for post-TCM-18 OD solutions. 

Over the course of the post-TCM-18 OD solutions, the nominal predicted landing location is observed to be sta- 
ble with little drift. The nominal predicted landing locations are hard to differentiate, as they lie nearly on top of 
each other. The benefit of stable OD solutions is that greater confidence can be placed on the footprint size refine- 
ment. For the post-TCM-18 OD solutions, the 99 percentile footprint ellipses decrease in size. The footprint size 
decreases from 113.3 km in downrange for OD s06011d to 75.5 km for OD s06014a. Table 3 summaries the varia- 
tion in the 99 percentile landed footprint ellipse sizes. The corresponding inertial flight-path angle errors at atmos- 
pheric interface for these three OD solutions show a similar trend and are also summarized in Table 3. 


7 

American Institute of Aeronautics and Astronautics 


Table 3: Variation in Landing Footprint Ellipse Size 


OD Solution 

99 Percentile 
Footprint Size (km) 

3-o Flight-Path 
Angle Error (deg) 

Post-TCM-18 

s0601 Id 

113.3 x 19.7 

±0.144 

s06012b2 

93.2 x 19.5 

±0.115 

s06014a 

75.5 x 19.4 

±0.083 

Post-TCM-19 

s06014c 

52.7 x 19.3 

±0.033 

s06015a 

47.5 x 19.2 

±0.021 

s06015f 

45.1 x 19.2 

±0.0017 

s06015f with E-2 hr balloon data 

38.1 x 8.2 

±0.0017 


Similarly, Monte Carlo analyses were performed for all the post-TCM-19 OD solutions (OD s06014b through 
OD s06014f). Figure 9 depicts the landing locations for a few of these OD solutions (OD s06014c, OD s06015a, and 
OD s06015f). All the OD solutions post-TCM-19 were extremely stable and produced nearly identical landing loca- 
tions as observed in Fig. 9, where the 99 percentile landing ellipses lie nearly on top of each other. The results for 
OD s06015f (which was the last OD solution available prior to entry) showed that all the EDL requirements and 
constraints were well within their limits, and that the final predicted nominal landing location was very close to the 
desired target (only 5.7 km away). The footprint size for these OD solutions decreases from 52.7 km in downrange 
for OD s06014c to 45.1 km in downrange for OD s06015f. Table 3 summaries the landed footprint variation, along 
with the variation in the corresponding inertial flight-path angle error. As seen, the inertial flight-path angle for OD 
s06015f was extremely small having a value of ±0.0017 deg. As such, the Navigation Team accurately delivered the 
capsule to the desired entry conditions. Based on OD s06015f results of the predicted nominal landing location be- 
ing 5.7 km away from the desired target and with a 99 percentile footprint ellipse of 45.1 km by 19.2 km, the 
authorization for capsule separation and subsequent Earth entry was granted. 
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Figure 9. Landing locations for post-TCM-19 OD solutions. 

Unfortunately, on entry day of January 15, 2006, a winter storm was moving through western Utah, which was 
producing very strong winds over UTTR. The effect of strong winds would cause the capsule during parachute de- 
scent to drift from its predicted landing location. To address such a scenario, the Stardust Entry Operations strategy 
had baselined two balloon launches (at E-8 hours and E-2 hours) to obtain measurements of the actual winds that 
would be observed during the entry in an effort to better predict the landing location and footprint ellipse. These two 
balloon measurement data revealed that very strong sustained winds were present over UTTR due to the winter bliz- 
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zard. Figures 6 and 7 show the actual wind profiles measured at E-2 hours. Note, although the balloon measurement 
data had a timestamp of E-2 hours, the actual measurements were taken approximately at E-4 hours due to the logis- 
tics and processing time requirements. As seen in Fig. 6, a very strong sustained wind to the North was present at E- 
2 hours. This sustained Northward wind component had a peak value of 50 m/s, and corresponded to approximately 
a 3 -a high profile from the GRAM-95 model. Although, the storm winds were predicted to subside over the remain- 
ing two hours prior to landing, this sustained Northward wind would tend to push the capsule landing location to- 
wards the North. Consequently, an updated prediction of the nominal landing location (using this balloon measured 
wind data) was necessary to aid the retrieval of the capsule by notifying the recovery team of the change in nominal 
landing location. 

Using this E-2 hour balloon measurement wind data, the Monte Carlo dispersion analysis was repeated with the 
GRAM-95 wind dispersions being replaced with the E-2 hour wind profile. Figure 10 shows the updated predicted 
nominal landing location and 99 percentile footprint ellipse. As seen, the OD s06015f landing prediction shifted 
North due to the E-2 hour balloon measurement wind data. The predicted nominal landing location moved from 5.7 
km West- Southwest of the target (using the GRAM-95 wind profiles) to 9.4 km North-Northwest of the target. The 
updated 99 percentile footprint ellipse was 38.1 km by 8.2 km. The Northward wind shifted the predicted nominal 
landing location approximately 1 1.5 km due North. The actual final landing location is shown in Fig. 10 and was 8.1 
km North-Northwest of the desired target, and was within the OD s06015f pre-entry predicted 99 percentile landing 
ellipse. Although, the actual final landing location indicated that the capsule had not drifted as much to the North as 
the updated prediction, this outcome was undoubtedly due to the winds subsiding (as forecasted) from that balloon 
measured data at E-2 hour, which was used in the update prediction. 

Target Location: (246.55 deg, 40.3167 deg) 



East Longitude, deg 

Figure 10. Final capsule landing location. 

VI. Conclusion 

On the morning of January 15, 2006, the Stardust entry capsule successfully descended through the Earth’s at- 
mosphere decelerating with the aid of a parachute and landed at the Utah Test and Training Range completing a 
7 year journey that returned cometary samples from the comet Wild-2. The navigation and the entry, descent, and 
landing trajectory analyses that were performed during the Mission Operations Phase upon final approach to Earth 
accurately delivered the entry capsule to the desired landing site. The capsule landed 8.1 km from the desired target, 
and was well within the allowable landing area at the Utah Test and Training Range. As a result, the process of tar- 
geting a capsule from deep space and accurately landing it on Earth was successfully demonstrated. Preliminary 
indications are that the actual entry trajectory was very close to the pre-entry prediction. As a result, the design prin- 
ciples and methodologies utilized for the entry flight dynamics, aerodynamics, and aerothermodynamics analyses 
were corroborated. 
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